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Gravitational Radiation from Triple Star Systems 

Enrico Montanari*, Mirco Calura, and Pierluigi Fortini 
Department of Physics, University of Ferrara and INFN Sezione di Ferrara, Via Paradiso 12, 1-44100 Ferrara, Italy 

We have studied the main features of the gravitational radiation generated by an astrophysical 
system constituted of three compact objects attracting one another (only via gravitational interac- 
tion) in such a manner that stable orbits do exist. We have limited our analysis to systems that can 
be treated with perturbative methods. We show the profile of the gravitational waves emitted by 
such systems. These results can be useful within the framework of the new gravitational astronomy 
which will be made feasible by means of the new generation of gravitational detectors such as LISA 
in a no longer far future. 
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I. INTRODUCTION 



A great deal of efforts are aimed at the detection of gravitational waves. Several large research programmes are 
under development all over the world flUEf . Consequently the theoretical study of gravitational radiation generated 
by astrophysical systems is quite a relevant topic. So far this study has been mainly focused on such sources as 

fvj ■ supernovae explosions, pulsars and compact binaries. In particular these latter sources have been widely investigated 
throughout all the stages of their evolution [||-[7| ■ 

Space based interferometric detectors such as LISA are sensitive to gravitational waves in the frequency band 1CP 4 - 
1CP 1 Hz j8|,|9| . Many known gravitational waves sources, such as short-period binary star systems, are contained in 
this bandwidth ll^JIlt]. However their waveforms could be biased by spurious effects such as those due to tidal forces, 
spinning of the stars, post-Newtonian corrections, presence of interstellar medium or of a third body. In order to 
recognize a gravitational wave from the experimental data, the expected wave form of the signal must be known 

f^ ■ in a highly accurate way. It is for this reason that many results have been obtained in order to take into account 
some of the above effects |L^-|l7|l. However, to our knowledge, a solution to the problem of the waveform from 

CT^ . triple star system is still lacking. In this paper we study the characteristics of the gravitational radiation produced 
by such a source, emphasizing the substantial differences with the waves originating from a simple binary system, 
^j*. providing a tool in order to recognize astrophysical systems from the analysis of the gravitational signal detected. 
Future gravitational astronomy will allow to open a new window on the universe besides the usual way to investigate 
(^Tj, the sky with electromagnetic devices or the relatively new neutrino astronomy (whose birth could be the "neutrino 
J* ' description" of the supernova explosion 1987A |lq-|2l|l). In fact gravitational radiation detection could be the only 
way to discover and study all that astrophysical systems which are not visible with standard methods because they 
do not emit electromagnetic radiation nor neutrinos. In particular, the existence of a black-hole could be argued by 
its effect on the gravitational radiation emitted by a binary system. In this paper we show that the gravitational 
signal of a binary system is modulated by a perturbing third body. We also show that the signal is different in form. 
Besides, triple star systems formed by a close binary and a third body radiate the same amount of energy of a close 
binary (or even more) in the region of interest for LISA. This means that this detector can have the same sensitivity 
for both kinds of system provided that the waveform of the radiation from triple star systems is known. In fact, with 
this knowledge, one can analyze the two time series that will be the data of LISA, searching for periodic sources not 
having the typical form of the signal from a binary system |9| . This can be done without any change in the detector 
project, but only in the data analysis. In this way it would be possible the gravitational analogue of the discovery of 
triple star systems with electromagnetic timing J22J. 

II. THE EQUATIONS OF MOTION 

Let us consider a system composed of three bodies attracting one another via gravitational force and assume that 
their typical sizes are small compared with their mutual distances. Let t*j and m* be the radius vector and the mass 
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of the i-th body in a given Galilean reference frame. In order to simplify the problem the following set of variables is 
introduced: 

miri + m. l r2 + m r 3 

xu) =r 2 -r 1 ; x, \ - r 3 - n; av cm ) = ■ ■ 

v ' mi + rrii + m 

m^ = m 2 ; m = m 3 . (2.1) 

where x^x) (A = i, o) are the relative coordinates of the A-th body with respect to the first one, while ai( cm ) are the 
coordinates of the centre of mass. The equations therefore become (all through the paper G = c = 1): 

(2.2) 
(2.3) 

where the so called disturbing functions have been defined 

^t^,^, - «. ( 1 _l_ I -=^a) (2.5, 

5£ e (iE(i),a;( )) = — - 6Ci(x( o) ,X(i)) (2.6) 

TTIq 

We assume that the mutual interaction between mi and m is a perturbation to their Keplerian motion around mi. 
This is accomplished if the following conditions are fulfilled: 

|a?(*)| « \x(o)\ (2.7) 

mi+m '»M 3 (2.8) 

m \X(o)J 

-J*—«(?V) 7 (2.9) 

mi + TO \X( )J 

Within the framework of these assumptions it is possible to solve the equations of motion by employing the so 
called method of variation of parameters in the form due to Lagrange. A complete review of this method, both in its 
classical formulation and in the post-Newtonian extension, may be found in Refs. p3J-p5| . 

In celestial mechanics the usual way to write the solution as a Keplerian motion is p3j: 

x\ = (cos(oj\) cos(Qa) — cos(u) SUi(Oa) sin(wA)) x\ + 

— (sin(cjA) cos(^a) + cos(za) sin(fi.\) cos(wa)) x\ (2-10) 
x\ = (cos(wa) sin(riA) + cos(za) cos^a) sin(wA)) x\ + 

— (sin(cjA) sin(fiA) — cos(i\) cos(Oa) cos(ui\)) i\ (2-11) 
x\ = sin(iA) sin(wA) i\ + sin(?A) cos(wa) i\ (2.12) 



where 



x\ = a x (cos(?7a) - e A ) (2.13) 



i\ = a \\J 1 - e l sin(?7A) (2.14) 

M x = n x (t - T x ) = 7 lx - e x sin(r)x) (2.15) 



/mi+m A , . 

UA = y^3— (2-16) 

(A = i, 6). Solutions are parameterized by means of the 12 Keplerian elements ax, e x , M x , wx, Qx, i\ p3| , p4| . Within 
the fr amew o rk of the method of variation of the parameters, the solution to the full problem has the same form as 
Eqs. ( 2.10 )-( fTl6 ) but with elements slowly varying upon time according to Lagrangian planetary equations |23|]24| ]. 



III. THE GENERATED GRAVITATIONAL WAVES 

The formal solution-we have reviewed in the previous Section-allows us to achieve the expression of the emitted 
gravitational wave as a function of the orbital elements and time explicitly. Since orbital parameters are slowly 
varying, we neglect their time derivative in performing the calculation. First we evaluate the reduced quadrupole 
moment of the whole system. According to its usual definition M we have: 
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where, in our case, 
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Substituting eq. (2.1) into the above relations we get 
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a and /? run from 1 to 3. The first two terms in Eq. ( p.2|) are the quadrupole moments of the two binary systems 
formed by mi and rrii (the inner one) and by m\ and m (the outer one); it retains the usual form but orbital elements 
vary with time. The third term, -an interaction one-is peculiar of the system as a whole. To clarify these points let us 
consider the components of the quadrupole tensor when eqs. (2.10-2.14) are taken into account. In general one has 



^ = E<(^>^^)^(^a) 
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^M = H d°f(u>i,tti, k, uj , n o , i )Qf (xj,x?,x l , x 2 ) 
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where df are related to the Euler matrices of the rotation linking the intrinsic reference frames of motion to the 
observer reference frame; Q u are the "intrinsic" quadrupole moment coefficients. If the orbital elements did not 
change, the first expression would be the usual q uadr upole tensor of a binary system. The effect of a third body thus 
induces a modulation of this signal. As for Eq. (3.4), it is a peculiar feature of a three body problem; however the 
structure of such a term is similar to the preceeding ones. 

The knowledge of the quadrupole tensor allows us to write the TT-gauge components of the gravitational waves 
emitted by the triple system. As it is well known (e.g. P,|26[|): 



. TT 



2 

3W 



d 2 D, 
dt 2 



— run 



d 2 D 



ok 



1 



S i] n l n k 



d 2 D, 



Ik 



(It 2 



(It 2 



-riinjn n 



, d 2 D l} 
d 2 D lk 



(It 2 



(3.5) 



where n — r/r is the unit vector in the observer direction. 

In order to obtain analytic expressions for the metric perturbation, it is simpler to consider the so called true 
anomaly ip\ instead of the eccentric anomaly rj\. The relation between these two angles is given by fl23f: 
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Starting from the above relation and the definition of rj\ [see eq. (2.15)] it is easy to show that the time derivative of 
ip\ is given by 
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In order to put the result in a compact form we introduce the Euler Matrices A\ with which eqs. ( 2.10| - ]2Tl2| ) can 
3 written as (here ' means transposition) 

xx = A\x x (3.6) 

The second derivative of the quadrupole tensor components can be written as (H G {11, 12, 0}) 

■D(l) = 3^ A ^d| J A)N (w A ,ft A ,i A ) J ff A (a A ,e A ,'i/'A) 
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where 
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The above coefficients depend only upon the parameters ax, e\ and M\, while 

(3.16) 
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(3.18) 
depend only upon the Euler angles u>\, fix and i\. The TT-gauge components of the gravitational field given in 
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Eq. (3.5) can then be written as 

h% T) =h$ +f tf + h^ (3.20) 

where 
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The values of the coefficients cA^ux, Qx,ix, n) and c, 3 io *. rs (u>i, 0J o ,Q,i,Q, o ,ii,i o ,n), connecting the Euler angles with 
the observer-source direction, are given in appendix. 

The signal is therefore the sum of three terms: two binary-like and an interaction one. As for the first two 
terms, they are equal in form to the signal emitted by two binary systems M; however the orbital elements are not 



constants; this fact results in a modulation of the signal. The third term in Eq. (3.20) is a peculiarity arising from 



the interaction between these two binary systems. Therefore the signal is described by a Fourier series with slowly 



variable coefficients. Within the framework of conditions (2.8)— (pT9[) the interaction term is of the same order of or 



greater than the inner one. Therefore, when the outer signal can be neglected (this may happen, for instance, if the 
sensitivity of the detector near n is negligibly small; see following section) , the signal form is different from that one 
emitted by a simple two-body system. 

IV. APPLICATIONS 

The time variation of orbital elements, induced by the disturbing functions, is very slow and it is quite difficult 
to observe within an interval of time which is of the same order of magnitude of the smaller unperturbed period 
involved. Therefore we are interested in the so called "secular" variation of the parameters. For the sake of simplicity 
we set %i = i = i, that is to say we assume the two orbital planes to coincide. Consequently we have fl; = G = Q 
(because t here is onl y on e line of nodes) . As we are only interested in secular variations of the elements, the disturbing 



functions (2.5) and (2.6) may be averaged to obtain: 
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This result is exact in both e» and e D ; see Ref. [£4J for an expansion up to second order in both eccentricities. By 
using Lagrange equations ||23| , [24f , up to first order of approximation, one gets: 



U>i = v i t J r LUio 

LJ = V t + Ld o0 

Mi = N l t + o l 
M = N t + a 

where we have put u)j = LUi + £1;, uj = lo + Q . Wjq, & q, &i and o~ are constants which are determined 
initial conditions are given, while v\ are the rates of periastron precessions. Their values are given by the 
relations: 
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By the way, we point out that the averaged Lagrangian given in Ref. p4( ], p. 320, does not yield the correct secular 
variation up to second order in eccentricities. In fact, inspection of Lagrangian equations shows that the averaged 
disturbing function should be expanded up to second order in both eccentricities, without neglecting terms proportional 
toefe 2 . 

Summing up, the secular effects of the interaction between m, and m are a change in the frequency of the 
unperturbed motion and appearanc e of t wo n ew frequencies related to the periastron precession. 

In the framework of assumptions (2.8)-( f2~9| ), which are the conditions for the perturbative approach to be allowed, 
we consider two possible systems. The first one is a close binary system orbiting around a third body whose mass 
is much greater than the total mass of the other two bodies. In other words the order of magnitude of the binary 
system size is much smaller than the mean distance between its centre of mass and the third body (the so called 
"lunar case"). 

The second system is a body orbiting around a close binary system whose total mass is much greater than that of 
the far body (the so called "solar system case" ) . 



A. First case 



Let us consider two stars whose masses are mi and mi << mi that, besides a motion around their centre of mass, 
are orbiting around a third star with a mass m » mi + m^. This two Keplerian orbits, called respectively inner 
and outer orbit, are perturbed by the disturbing functions (4.1)— (4.2). The effect of these perturbative terms, as we 



have already seen, is to cause a variation of the elements of the unperturbed Keplerian motions. 

In order to show the peculiarity of the gravitational signal in such a situation we consider the following example 
in which we have assigned: mi = 2 • 10 31 Kg « 10 m©, m, = 2 • 10 29 Kg « 10 _1 m©, m = 2 • 10 32 Kg « 100 m©, 
Oj = 10 9 m, a — 10 Oj |H| = 10 kpc, ei = e e = 0.25, i = 0, CI = 0, n = (0, 0, 1), w i0 = 1.1, a) o0 = 0.4, Oi = 0.3 and 
(t = 0.1. It could be formed by a neutron star close to a black-hole a solar diameter apart, under the influence of a 
farther more massive black-hole. 

The motion of this system could be viewed, in zero-th order of approximation, as the Keplerian motion of the 
centre of mass of the two closest stars around the far body and the relative motion, Keplerian too, of the binary 
system. The perturbation, due to the interaction among the three bodies, is therefore described by the variation of 
the elements of these two Keplerian orbits. The effect results in the appearance of two new frequencies (the rates of 



peria stron p rece ssion, s ee E qs, (4.7) and (4.8)) in the gravitational signal emitted besides the two Keplerian ones (see 
Eqs Q2.16D , (pi, and (|4.10|)). In this case we find 



Ni = 1.14 X 10- 3 — ; 7V o = 1.21 x 10- 4 — ; 

s s 

^ = 9.19xl0- 6 — : ^ O = 9.64xl0- 10 — 

s s 

We notice that a similar effect is due to the 1PN perturbation of the inner orbit; in this case the 1PN frequency is 
given by |J: 

v\ PN = 3(Wl + m ;} n t = 5.56 x 10- 8 rad/s. (4.11) 

<H (1 - ef ) 

Since our treatment of the phenomena is a classical one, its validity holds for observational times smaller than 
27r/V 1PAr ~ Ay. Therefore, for this kind of systems, only frequencies greater than this one, that is to say i%i, n and 
Vi could be appreciate. 

In the previous section we have shown that the gravitational wave emitted is the sum of three terms: the first 
(second) one is due to the inner (outer) motion and involves two frequencies, that is to say N% (N ) and Vi (v ); the 
last one is an interaction term which depends on all the four frequencies. From what was said it is clear that for the 
outer term the modulation of signal cannot be observed, while, as far as the interaction term is concerned, only two 
frequencies yield an observable modulation. 

A gravitational detector presents a sensitivity which depends on the frequency. If we consider a detector for which 
the ratio between sensitivity near N and near Ni is negligible small (such an assumption is quite realistic for LISA, 
see for instance [pi) then the detectable wave profile is merely the sum of the inner and interaction terms even if 
the greater amplitude is associated with the outer motion (as it is in this case). The output results in a signal with 
frequency Ni whose amplitude is modulated by the two frequencies N and i>i (inner precession rate). 

In the figures below (figs. 1-3) we show the gravitational wave form for the inner term, the interaction one and 
their sum. Time is given in units of the period of the inner orbit. 

B. Second Case 

The second kind of triple-system we consider, involves a star of mass toi and other two bodies (whose masses are 
rrii and m ) moving around it; in addition the following relation holds: m\ 3> m 3> rrii. We assume that the motion 
of both rrii and m around the central body mi is a Keplerian one with orbital elements varying with time because 
of the mutual interaction between the two lighter bodies. 

Now, as it has already be done in the previous case, we study the gravitational waves generated by such a kind 
of system in a particular case. In doing so we have set: mi — 2 x 10 32 Kg « lO 2 m0, m, = 2x 10 29 Kg « 10 _1 to q , 
to = 2 x 10 31 Kg « lOm , a,=4x 10 9 to, a = 10 a t , \R\ = 10 kpc e t = e e = 0.25, i = 0, Q = 0, n = (0,0, 1), 
ujio = 1.1, (I> o = 0.4, <j, = 0.3 and a = 0.1. In this case the two Keplerian frequencies of the inner and outer orbit 
(Ni and N respectively) are: 

7V 4 = 4.57 x nr 4 — ; N = 1.52 x lO" 5 — (4.12) 

s s 

while the two periastron precession rates are: 

^=3.65xl0- 8 — ; v = 1.20 x lO" 10 — (4.13) 

s s 

The frequency involved in the 1PN perturbation of the inner orbit is: 

v]PN = 3(toi + to,) m 5 ^ x iQ _ s rad 
at (1 - ef) s 



From (4.14) it follows that the validity of our results holds for times smaller than 2-kJv\ pn = 1.2 x 10 8 sec ~ Ay, 
therefore we can not appreciate the frequencies Vi and v . For this reason the detectable frequencies of the signal are 
the only Keplerian ones; it follows that the inner and outer term of the wave profile are not modulated, while the 
interaction term is modulated by only one frequency. 

As in the first case, we show graphically (figs. 4-5) the gravitational wave form; the observation time is measured 
in units of the periods of the inner orbit. Also in this case if the ratio between the detector sensitivities near N and 
Ni is negligible small, then the detectable signal is the sum of the inner and interaction terms, which gives a signal 
with frequency Ni and amplitude modulated by the outer frequency N (fig. 6). 



V. CONCLUSIONS 

In this paper we have studied the gravitational radiation emitted by triple-star systems. 

We have shown that the signal emitted by a triple star system such as that envisaged here is the sum of three 
terms. Two of them are equal in form to the signals from two binary systems, each of them modulated by the secular 
variation of the orbital elements. The last term originates in the three-body interaction, not having a two-body 
counterpart. It has a fundamental frequency near m (angular frequency of inner motion) and it is modulated by n 
(angular frequency of outer motion) and by the secular variation of orbital elements. If rii is the only frequency in the 
range of detector sensitivity (for LISA such an assumption is quite realistic) then a close binary system in interaction 
with a third body, in the assumptions made through the paper, gives rise to an output signal which is different in 
form from the signal belonging to an isolated two-body system. 

If the inner motion is in the region of this detector sensitivity, the waveform could be experimentally extracted by the 
two time series which will be the data of the space based interferometer LISA. Therefore space based interferometers 
able to detect gravitational radiation from galactic close binary systems will make possible the discovery and study 
of eventual three-body systems, without any change in their planned design. 

In order to test the applicability of this approach we have considered two particular cases. We have found that, 
because of the Newtonian approximation we used, there could be systems of interest for which this treatment is enough 
precise to answer for three (first case) or two (last case) modulating frequencies. In this paper we have assumed the 
post-newtonian effects involved in motion of the bodies to be negligible. Actually, were two of the three bodies very 
close to each other, such effects would become too important to be neglected. In this case our classical perturbative 
approach to the motion should be replaced by its post-Newtonian extension (see p5j, where relativistic corrections 
in a simple three-body system are investigated). 
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APPENDIX: 



In the following we give the values of the coefficients c,°y-,^{u!\, tt\,i\, n) and cP. •, _(u>j, u , 0,i,Q o ,ii,i o , n) which 

enter in eqs. ( 3.21 )— ( 3.22 ). For simplicity we have omitted lower indices, being understood that the following relations 
hold for every value of them. 
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FIG. 1. Amplitude of the term h^ of the signal emitted by the triple star system o f Sec tio n 4. 1; the orbital plane is 



perpendicular to the observation direction. For the meaning of the symbols used see eqs. (3.20)- (3.22). a) The fundamental 



frequency is Ni/-K. b) The modulation of the signal has an angular frequency given by the rate of periastron precession 8m. 
FIG. 2. Amplitude of the interaction term n£° emitted by the triple star system o f Sec tio n 4.1 ; the orbital plane is 



perpendicular to the observation direction. For the meaning of the symbols used see eqs. ( 3.20 )- ( 3.22 ). a) The fundamental 
frequency is A^/(2 7r) while the frequency modulation is N /(2-k). b) The third modulating angular frequency is the periastron 
precession rate 5rn. 

FIG. 3. Amplitude of n{° + h^l, the detectable signal (see considerations in the text) emitted by the triple star system of 
Section 4.1. a) It is the sum of figs, la and 2a. b) It is the sum of figs, lb and 2b. 



FIG. 4. Amplitude of the inner term h{{ of the gravitational wave emitted by the triple star system of Se cti on 4 .2; the 
orbital plane is perpendicular to the observation direction. For the meaning of the symbols used see eqs. (3.20)- (3.22). The 
fundamental frequency is Ni/n. 



FIG. 5. Amplitude of the interaction term h}[° emitted by the triple star system of Secti on 4 .2, 
perpendicular to the observation direction. For the meaning of the symbols used see eqs. ( 3.2C )— 
frequency is Ni/(2ii)\ the modulation is given by N /(2n). 



when the orbital plane is 
J3.22|). The fundamental 



FIG 
Section 4.2 



6. Amplitude of n{° 



+ /ijj , the detectable signal (see considerations in the text) emitted by the triple star system of 



It is the sum of figs. 4 and 5. 
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